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Abstract. The Fat-Link Irrelevant Clover (FLIC) fermion action provides a new form 
of nonperturbative 0(a)-improvement in lattice fermion actions offering near contin- 
uum results at finite lattice spacing. It provides computationally inexpensive access 
to the light quark mass regime of QCD where chiral nonanalytic behaviour associated 
with Goldstone bosons is revealed. The motivation and formulation of FLIC fermions, 
its excellent scaling properties and its low-lying hadron mass phenomenology are pre- 
sented. 

1 Introduction 

The origin of the masses of light hadrons represents one of the most fundamental 
challenges to QCD. Despite the universal acceptance of QCD as the basis from 
which to derive hadronic properties, there has been slow progress in understand- 
ing the generation of hadron mass from first principles. Solving the problem of 
the hadronic mass spectrum would allow considerable improvement in our un- 
derstanding of the nonperturbative nature of QCD. The only available method 
at present to derive hadron masses directly from QCD is a numerical calculation 
on the lattice. 

The high computational cost required to perform accurate lattice calculations 
at small lattice spacings has led to increased interest in quark action improve- 
ment. In this article we present results of simulations of the spectrum of light 
mesons and baryons using an 0(a) improved fermion action [1-4]. In particular, 
we will start with the standard clover action and replace the links in the irrele- 
vant operators with APE smeared [5,6], or fat links. We shall refer to this action 
as the Fat-Link Irrelevant Clover (FLIC) action. Although the idea of using fat 
links only in the irrelevant operators of the fermion action was developed here 
independently, suggestions have appeared previously [7] . 

In Section 2, we provide the reader with some background information on lat- 
tice fermion actions. In particular, we start with the basic lattice discretisation 
of the derivative appearing in the continuum Dirac action, followed by improve- 
ments suggested first by Wilson [8] and then by Sheikholeslami and Wolhert 
[9] . Section 3 contains the procedure for creating the FLIC fermion action while 
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in Section 4 we describe the gauge configurations used in our lattice simula- 
tions. The results of an investigation of the scaling of this action at finite lattice 
spacing are presented in Section 5. In Section 6 we investigate the problem of 
exceptional configurations by performing simulations of hadron masses at light 
quark masses corresponding to mT,/mp = 0.35. Section 7 discusses the evidence 
for enhancement in octet-decuplet mass splittings as one approaches the chiral 
limit in the quenched approximation and finally in Section 8 we summarise the 
results. 



2 The Lattice Quark Action 
2.1 The Naive Fermion Action 

In Euclidean space-time, the continuum Dirac action is written as 

v;(^+m)v, (1) 

where the covariant derivative is defined as = dfj, + i g A^. Wilson [8] discre- 
tised the continuum Dirac action by replacing the derivative with a symmetrised 
finite difference. Gauge links are included to not only encode the gluon field, 
but to also maintain gauge invar iance 



Ufj,{x) ip{x + ajl) — UUx — afb) tp{x — ajl) 



(2) 



Our notation uses the Pauli (Sakurai) representation of the Dirac 7-matrices de- 
fined in Appendix B of Sakurai [10] . In particular, the 7-matrices are hermitian, 

7p = 7^ and cr^,^ = [7^, ■y^]/{2i) such that 7^7^ = +icr^,.. The gauge link 
variables at space-time position x are defined as 



Ui^{x) = V exp g J A^(a; -h A/i) rfAj 



(3) 



Here the operator P path-orders the A^^'s along the integration path, a is the 
lattice spacing, and g is the coupling constant. 

The continuum Dirac action is recovered in the limit a ^ by Taylor ex- 
panding the t/^ and ip{x+afi) in powers of the lattice spacing a. Keeping only the 
leading term in a (and for ease of notation we write a/i — > /t), Eq. (2) becomes 

-^V^(a;)7^[ ( l + iagAi,{x+^) + ...){i>{x)+a^\x) + ...)- 

( 1 -iagA^ix - |) + • • ■){i^{x) - a^'(a;) + ...)] 
= Vi(a;)7^ + 0{a''))i;{x)+igi;ix)j^[A^ + 0{a^)]i;{x), (4) 



which is the kinetic part of the standard continuum Dirac action in Euclidean 
space-time to O(a^). Hence we arrive at the simplest ("naive") lattice fermion 
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action, 

Sn = mg^ tp{x) tp{x) + 

X 

u^{x) ip{x + A) - ul{x - A) i}{x - fi) 

= ^^(x)Mi;[C7]V(y), (5) 

X 

where the interaction matrix is 

The Taylor expansion in Eq. (4) shows that the naive fermion action of Eq. (5) 
has 0{a^) errors. It also preserves chiral symmetry. However, in the continuum 
limit it gives rise to 2^^ = 16 flavours of quark rather than one. This is the famous 
doubling problem and is easily demonstrated by considering the inverse of the 
free field propagator (obtained by taking the fourier transform of the action with 
all U^ = l) 

S~'^{p) = mq + -'Y^^^smp^a, (7) 

which has 16 zeros within the Brillouin cell in the limit mq 0. eg, = 
(0,0,0,0), (tt/o, 0,0,0), (7r/a,7r/a, 0, 0), etc. Consequently, this action is phe- 

nomenologically not acceptable. 

There are two approaches commonly used to remove these doublers. The 
first involves adding operators to the quark action which scale with the lattice 
spacing and thus vanish in the continuum limit. These operators are chosen to 
drive the doublers to high energies and hence are suppressed. This technique 
for improving fermion actions proceeds via the improvement scheme proposed 
by Symanzik [11] and is discussed in more detail in the following sections. The 
second method for removing doublers involves "staggering" the quark degrees of 
freedom on the lattice. This procedure exploits the fact that the naive fermion 
action has a much larger symmetry group, C/y(4)^ f/4(4), to reduce the doubling 
problem from 2'' = 16 ^ 16/4 while maintaining a remnant chiral symmetry. 
This approach is not used in this article so the details of the action will not be 
discussed here. Details of the derivation of staggered fermions can be found in 
most texts (eg. [12,13]). 

2.2 Wilson Fermions 

In order to avoid the doubling problem, Wilson [8] originally introduced an 

irrelevant (energy) dimension- five operator (the "Wilson term" ) to the standard 
naive lattice fermion action (Eq. (5)), which explicitly breaks chiral symmetry 
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at 0{a). 



2 ^^^^^ 



m 



where r is the "Wilson coefficient," 

V^^(a;) = [C/^ (a;) ij{x + fi) - U^ix - fi) ^p{x - jl)] , 

and 



(8) 



(9) 



= ^[U^i^) + A) + Ulix - A) V(a; - A) - 2 (10) 

The Wilson action in Eq. (8) has no doublers for r > as the Wilson term gives 
the extra fifteen species at = n/a a mass proportional to r/a. Also, if r = 1 
then the Wilson action has no ghost branches in its dispersion relation (see, for 
example, Ref. [14]). In terms of link variables, U^{x), the Wilson action can be 
written 



Too 



4r \ - 1 - 



-{in + t) Ul{x - A) i){x - jl) 
where the interaction matrix for the Wilson action, M^, is usually written 



M. 



w 



(11) 
(12) 

1 

(13) 



with a field renormalisation 

K = 1/(2to5 a + 8r) , 
= V/a/2k. 

We take the standard value r = 1 and the quark mass is given by 

1/1 1 



(14) 



(15) 



In the free theory the critical value of kappa, Kc, where the quark mass vanishes, 
is l/8r. In the interacting theory. Kc is renormaliscd away from l/8r. The quark 
mass has both multiplicative and additive renormalisations due to the explicit 
breaking of chiral symmetry by the Wilson term. 
In the continuum limit, one finds 

Sw = Jd^xxPix) (^]/)+m-^Y.Dlj t/-(x)+0(a). (16) 
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By lifting the mass of the unwanted doublers with a second derivative, 0{a) 
discretisation errors have been introduced into the fcrniion matrix. In contrast, 
the Wilson gauge action [8] has only 0{a'^) errors. Hence there is enormous 
interest in applying Symanzik's improvement program [11] to the fermion action 
to remove 0{a) errors by adding additional higher dimension terms. 

2.3 Improving The Fermion Action 

The addition of the Wilson term to the fermion action introduces large 0{a) 
errors which means that in order to extrapolate reliably to the continuum limit, 
simiilations must be performed on fine lattices, which arc therefore very com- 
putationally expensive. The scaling properties of the Wilson action at finite a 
can be improved by introducing any number of irrelevant operators of increasing 
dimension whose contributions vanish in the continuum limit. 

The first attempt at removing these 0{a) errors was by Hamber and Wu [15] 
who added a two link term to the Wilson action 



The removal of the classical 0{a) terms is easily observed through a Taylor ex- 
pansion. While this action also removes 0{a^) errors at tree-level, it has only 
received a small amount of interest due to the computational expense in evaluat- 
ing the double hopping term. Calculations that have been done with this action 
show that it works well at coarse lattice spacings and has the added bonus that 
it has an improved dispersion relation [16]. 

A more systematic approach [11] to 0(a) improvement of the lattice fermion 
action in general [9] is to consider all possible gauge invariant, local dimension- 
five operators, respecting the symmetries of QCD 



Shw=Sw + ^^(^) 




l^' + l^i') ^''(^) ^/'(^ + ^''(^ + 2A) 




i g a Csw r 
4 




O3 



bganiq 
2 



(18) 




Note that the operator Ipp is linearly related to Oi and O2 as 



p]l)=D^D^-^- 




(19) 
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Operator Oi is a new local operator in the lattice fermion action and must be 

included. On the other hand, O3 and Os of Eq. (18) act to simply renormalise 
the coefficients of existing terms in the lattice action, removing 0{amq) terms 
from the relation between bare and renormalised quantities. For example, the 
renormalisation of the quark mass rUq — > mq{\ — bmamq) incorporates O5. 
Similarly, O3 introduces a mass dependence in the gauge coupling ^ {1 — 
bganiq) such that the lattice spacing remains constant for constant g as rUq is 
varied [17]. For the quenched approximation, bg = 0. 

The key observation to efficient 0{a) improvement is that the 0{a) improve- 
ment afforded by two-link terms of the fermion action [18] may be incorporated 
to 0{a) into the standard Wilson fermion action complemented by Oi though 
the following transformation of the fermion fields 

ip ^ ip' = {1 + bqV a niq) {1 - Cqrap)^p , 

^p^Tp' = {l + bqramq)^{l + cqra'p), (20) 

where "0' represents the physical fermion field recovered in the continuum limit, 
while ip is the lattice fermion field used in the numerical simulations. The Jaco- 
bian of the transformation is 1 to 0{a) [9]. At tree- level, 6g = Cg = 1/4 correctly 
incorporates the 0{a) corrections of O2 and O4 into the fermion action. Note 
that this field transformation renormalises the fermion mass of the simulations 

m — > m(l -h ^ r am) . (21) 

For the spectral quantities investigated herein, it is sufficient to work with 
the lattice fermion-field operators alone. Here the fermion operators, act only 
as interpolators between the QCD vacuum and the state of interest and do not 
affect the eigenstates of the QCD Hamiltonian. However, for matrix elements 
of fermion operators, hadron decay constants, or off-shell quantities such as the 
quark propagator, it is important to take this field redefinition into account to 
connect the lattice field operators to the continuum field operators incorporating 
important 0{a) contributions. 

In summary, Oi, the "clover" term, is the only dimension- five operator ex- 
plicitly required to complement the Wilson action to obtain 0{a) improvement. 
This particular action is known as the Sheikholeslami-Wohlert fermion [9] action 

Ssw = Sw -'^-^^^^^^^ '>p{x)(Tf,i, Ff,^Tp{x) , (22) 

where Csw is the clover coefficient which can be tuned to remove 0(a) artifacts 
to all orders in the gauge coupling constant g. 

^ ^ f 1 at tree-level , 

\ I/uq mean- field improved , ^ ' 



with Uo the tadpole improvement factor correcting for the quantum renormali- 
sation of the operators (see definition in Section 4). Nonperturbative (NP) 0{a) 
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Fig. 1. Loops required to construct 



improvement [19] uses the axial Ward identity to tune Csw remove all 0{a) 
artifacts provided one simultaneously improves the coupling g^, the quark mass 
iriq, and the currents [19]. The advantage of the clover action is that it is local 
and is only a small overhead on Wilson fermion simulations. Further details of 
the improvement the clover action provides at finite lattice spacing is given in 
Section 5. 

The name "clover" is associated with the SW fermion action due to the lattice 
discretisation of the field strength tensor, i^^^. An expression for F^^, is obtained 
by considering the sum of the four 1x1 link paths surrounding any lattice site 
in the fi — v plane as shown in Fig. 1. Using the expansion for the elementary 
link product (see for example Ref. [12,13,20] ), we obtain the lattice expression 
for Ffj,^ 



(0W(:r)-0Wt(x) 



-^Tr(0W(x)-0«t(.) 



(24) 



where F^^ is made traceless by subtracting 1/3 of the trace from each diagonal 
element and 

O^^Kx) = U^{x) U,{x + fi) Ulix + 0) Ulix) 
+ U^{x) ulix + v- jj)Ul{x^ ^) Uf,{x - (l) 
+ Ul{x - fi) ulix ~ fi- v) Uf_,{x - /t - i>) U^{x - v) 
+ Ul{x - v) U^{x - v) U,{x + fi-v) ulix) . (25) 

Substantial progress has been made to improve -F^^ to 0{a^) by adding terms 
constructed using larger loops [20]. 
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3 Fat-Link Irrelevant Fermion Action 

The established approach to nonperturbative (NP) improvement [19] tunes the 
coefficient of the clover operator to all powers in . Unfortunately, this for- 
mulation of the clover action is susceptible to the problem of exceptional con- 
figurations as the quark mass becomes small. Chiral symmetry breaking in the 
clover fermion action introduces an additive mass renormalisation into the Dirac 
operator that can give rise to singularities in quark propagators at small quark 
masses. In practice, this prevents the simulation of small quark masses and the 
use of coarse lattices (/3 < 5.7 ^ a > 0.18 fm) [21,22]. Furthermore, the plaque- 
tte version of i^^y, which is commonly used in Eq. (22), has large 0{a'^) errors, 
which can lead to errors of the order of 10% in the topological charge even on 
very smooth configurations [23]. 

The idea of using fat links in fermion actions was first explored by the MIT 
group [24] and more recently has been studied by DeGrand et al. [21,25,26], who 
showed that the exceptional configuration problem can be overcome by using 
a fat-link (FL) clover action. Moreover, the renormalisation of the coefficients 
of action improvement terms is small. In principle it is acceptable to smear 
the links of the relevant operators. The symmetry of the APE smearing process 
ensures that effects are 0{a?). The factors multiplying the link and staple ensure 
the leading order term is e"'^'*'', an element of SU(3). Issues of projecting the 
smeared links to SU(3) are 0{a^) effects and therefore correspond to irrelevant 
operators [27]. However, the net effect of APE smearing the links of the relevant 
operators is to remove gluon interactions at the scale of the cutoff. While this has 
some tremendous benefits, the short-distance quark interactions are lost. As a 
result decay constants and vector-pseudoscalar mass splittings of heavy mesons, 
which are sensitive to the wave function at the origin, are suppressed [28]. 

The solution to this is to work with two sets of links in the fermion action. 
In the relevant dimension- four operators, one works with the untouched links 
generated via Monte Carlo methods, while the smeared fat links are introduced 
only in the higher dimension irrelevant operators. The effect this has on de- 
cay constants and vector-pseudoscalar mass splittings of heavy mesons is under 
investigation and will be reported elsewhere. 

Fat links [21,25] are created by averaging or smearing links on the lattice 
with their nearest neighbours in a gauge covariant manner (APE smearing). 
The smearing procedure [5,6] replaces a link, Ufj,{x), with a sum of the link and 
a times its staples 



4 




+ t/ 1 (,x- - I)) (7^ (x - z>) [/, (x - + A) J , (26) 

followed by projection back to SU(3). We select the unitary matrix U^^ which 
maximises 



netv{U^^U'J), 
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by iterating over the three diagonal SU(2) subgroups of SU(3). Performing eight 
iterations over these subgroups gives gauge invariance up to seven significant 
figures. We repeat the combined procedure of smearing and projection n times. 
We create our fat hnks by setting a = 0.7 and comparing n = 4 and 12 smearing 
sweeps. The mean-field improved FLIC action now becomes 



Si^ = S^^ - 'l^Z,Y ^(^) V'(x) , (27) 



where Ff^,^ is constructed using fat hnks, v^^ is calculated in an analogous way to 
Eq. (33) with fat links, and where the mean-field improved Fat-Link Irrelevant 
Wilson action is 



Uo Uo 



FL ^(^ + w + — — FL — '^y^ - w 



(28) 



As reported in Table 1, the mean- field improvement parameter for the fat 
links is very close to 1. Hence, the mean- field improved coefficient for Csw is 
accurate.^ It is in this way that the extensive task of non-perturbatively calcu- 
lating the renormalisations of the improvement coefficients discussed in Sec. 2.3 
is avoided. APE smearing the links of dimension five operators suppresses the 
renormalisation, allowing the precise matching of improvement coefficients with 
only tree-level knowledge of their values. 

In addition, one can now use highly improved definitions of F^^i, (involving 
terms up to mJ^), which give impressive near-integer results for the topological 
charge [20]. In particular, we employ the 3-loop C'(a'')-improved definition of 
F^v in which the standard clover-sum of four 1x1 loops lying in the /z, v plane 
is combined with 2x2 and 3x3 loop clovers. Bilson-Thompson et al. [20] find 



9 — Q 



(29) 

Tracclcss 



where W"^" is the clover-sum of four n x n loops and F^^ is made traceless 
by subtracting 1/3 of the trace from each diagonal element of the 3x3 colour 
matrix. This definition reproduces the continuum limit with 0{a^) errors. On 
approximately self-dual configurations, this operator produces integer topolog- 
ical charge to better than 4 parts in 10*. We also consider a 5-loop improved 
F^i, for the 20^ x 40 lattice at /? = 4.53. Since the results for the 5-loop operator 
agree with the 3-loop version to better than 4 parts in 10'' [20], we are effectively 
using the same action as far the scaling analysis is concerned. 

The use of thin links in Eq. (28) ensures that the relevant dimension-four 
operators sec all the dynamics of the Monte-Carlo generated gauge fields. Upon 

^ Our experience with topological charge operators suggests that it is advantageous to 
include wo factors, even as they approach 1. 
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Table 1. The value of the mean link for different numbers of APE smearing sweeps, 
n, at a = 0.7 on a 16^ x 32 lattice at /? = 4.60 which corresponds to a lattice spacing 
of 0.122(2) fm set by the string tension. 



n 


„ FL 
Mo 







0.88894473 


0.62445197 


4 


0.99658530 


0.98641100 


12 


0.99927343 


0.99709689 



expanding the thin links in terms of the gauge potential, 0{a^) contributions 
of energy- dimension six are revealed, which, ideally, should be removed via the 
fat-link irrelevant operator procedure. Fortunately, actions with many irrelevant 
operators {e.g. the D234 action) can now be handled with confidence as tree-level 
knowledge of the improvement coefficients is sufficient. However, as we will see, 
the scaling of the 0(a)-improved FLIC fermion action of Eqs. (27) and (28) is 
excellent already, and the added computational expense of the two-link hopping 
terms required in 0(a^)-improvement is not well motivated at present. 

Work by DeForcrand et al. [29] suggests that 7 cooling sweeps are required 
to approach topological charge within 1% of integer value. This is approximately 
16 APE smearing sweeps at a = 0.7 [30]. However, achieving integer topologi- 
cal charge is not necessary for the purposes of studying hadron masses, as has 
been well established. To reach integer topological charge, even with improved 
definitions of the topological charge operator, requires significant smoothing and 
associated loss of short-distance information. Instead, we regard this as an upper 
limit on the number of smearing sweeps. 

Using unimproved gauge fields and an unimproved topological charge opera- 
tor. Bonnet et al. [23] found that the topological charge settles down after about 
10 sweeps of APE smearing at a = 0.7. Consequently, we create fat links with 
APE smearing parameters n = 12 and a = 0.7. This corresponds to ^ 2.5 times 
the smearing used in Refs. [21,25]. Further investigation reveals that improved 
gauge fields with a small lattice spacing (o = 0.122 fm) are smooth after only 
4 sweeps. Hence, we perform calculations with 4 sweeps of smearing at a = 0.7 
and consider n = 12 as a second reference. Table 1 lists the values of Uq^ for 
n = 0, 4 and 12 smearing sweeps. 

We also compare our results with the standard Mean-Field Improved Clover 
(MFIC) action. We mean-field improve as defined in Eqs. (27) and (28) but with 
thin links throughout. For this action, the standard 1-loop definition of F^^, is 
used. 

4 Lattice Simulations 

The simulations arc performed using the Luscher-Weisz [31] mean-field improved, 
plaquette plus rectangle, gauge action on 12^ x 24 and 16^ x 32 lattices with 
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lattice spacings of 0.093, 0.122 and 0.165 fm determined from the string tension 
with = 440 MeV. We define 

= y E - f^«<i(^)) - ll^ E - f^rect(a;)) , (30) 

sq ^ rect 

where the operators Usqix) and Urecti^^) defined as 

f/s, (x) = U, (x) U„ {x + fi) Ul {x + v) Ul {x) , (31) 

Urect{x) = U^{x) Uy{x + jl)Uy{x + V + fl) 

xUl{x + 2i))Ul{x + v)Ul{x) 
+ Uf,{x) Uf,{x + (i) Uu{x + 2/t) 

X Ul{x + il + v) Ulix + y) Ul{x) . (32) 

The hnk product Urect{x) denotes the rectangular 1x2 and 2x1 plaquettes, 
and for the tadpole improvement factor we employ the plaquette measure 

/I 

uo= l-neti{U,q)\ . (33) 

Initial studies of FLIC, mean- field improved clover and Wilson quark actions 
were made using 50 configurations. The scaling analysis of FLIC fermions was 
performed with a total of 200 configurations at each lattice spacing and volume. 
In addition, for the light quark simulations, 94 configurations arc used on a 
20^ X 40 lattice with a = 0.134(2) fm. Gauge configurations are generated using 
the Cabibbo-Marinari pseudo-heat-bath algorithm with three diagonal SU(2) 
subgroups looped over twice. Simulations are performed using a parallel algo- 
rithm with appropriate link partitioning [32] , and the error analysis is performed 
by a third-order, single-elimination jackknife, with the per degree of freedom 
(A^df) obtained via covariance matrix fits. 

A fixed boundary condition is used for the fermions by setting 

Ut{x,nt)^Q and Uf^{x.nt) = {) V £C (34) 

in the hopping terms of the fermion action. The fermion source is centered at the 
space-time location [x, y,z,t) = (1,1,1,3), which allows for two steps backward 
in time without loss of signal, for all simulations except those on the 20"^ x 
40 lattice at /3 = 4.53 which has the fermion source located at {x,y,z,t) = 
(1,1,1,8). Gauge- invariant Gaussian smearing [33] in the spatial dimensions is 
applied at the source to increase the overlap of the interpolating operators with 
the ground states. The source-smearing technique [33] starts with a point source, 

iJoaix,t) = S'^'S^^S^^^Jt^to (35) 

for source colour c, Dirac 7, position xq = (1, 1, 1) and time to and proceeds via 
the iterative scheme, 

'^i{x,t) =^F{x,x')ipi-i{x',t) , 

x' 
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where 

Repeating the procedure TV times gives the resulting fermion source 

7PN{x,t)=J2F''ix^x')Mx',t). (36) 

x' 

The parameters N and a govern the size and shape of the smearing function. 
We simulate with N = 20 and a = 6, except on the 20^ x 40 lattice which has 
A'' = 35. The propagator, S, is obtained from the smeared source by solving 

M^'^Sf^ = ra, (37) 

for each colour, Dirac source c, 7 respectively of Eq. (35) via the BiStabilised 
Conjugate Gradient algorithm [34]. 



5 Scaling of FLIC Fermions 

Hadron masses arc extracted from the Euclidean time dependence of the calcu- 
lated two-point correlation functions. For baryons the correlation functions are 
given by 

Git;p,r) = ^e-'P- r^'^{n\T[x"ix)x''{0)m , (38) 

X 

where x are standard baryon interpolating fields, 12 represents the QCD vacuum, 
r is a 4 X 4 matrix in Dirac space, and a, (3 are Dirac indices. At large Euclidean 
times one has 

G{t; P,n^^^ e-^^' tr [Fi-ij -p + M)] , (39) 

where Z represents the coupling strength of x(0) to the baryon, and Ep = 
(p^ + M^)^/^ is the energy. Selecting p = and F = {1 + 74)/4, the effective 
baryon mass is then given by 

M{t) = log[G(t)] - log[G{t + 1)] . (40) 

Meson masses are determined via analogous standard procedures. The critical 
value of determined by linearly extrapolating as a function of iriq 

to zero. 

Figure 2 shows the nucleon effective mass plot for the FLIC action on a 
IG"* X 32 lattice at /? = 4.60 which corresponds to a lattice spacing of 0.122(2) fm 
set by the string tension. The fat links are created with 4 APE smearing sweeps 
at a = 0.7 ("FLIC4"). The effective mass plots for the other hadrons are similar, 
and all display acceptable plateau behavior. Good values of x^/^df are obtained 
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Fig. 2. Effective mass plot for tlic nucleoli for the FLIC action from 200 configurations 
on a 16^ X 32 lattice at /3 = 4.60 which corresponds to a lattice spacing of 0.122(2) fm 
set by the string tension. The fat links are created with 4 sweeps of smearing at a = 0.7. 
The five sets of points correspond to the k values listed in Table 2, with k increasing 
from top down. 



for many different time-fitting intervals as long as one fits after time slice 8. All 
fits for this action are therefore performed on time slices 9 through 14. For 
the Wilson action and the FLIC action with n = 12 ("FL1C12"), the effective 
mass plots look similar to Fig. 2 and display good plateau behavior. The fitting 
regimes used for these actions are 9-13 and 9-14, respectively. 

The values of k used in the simulations for all quark actions arc given in 
Tabic 3. We have also provided the values of Kci for these fcrmion actions when 
using our mean-field improved, plaquette plus rectangle, gauge action at /? = 
4.60. We have mean-field improved our fermion actions so we expect the values 
for to be close to the tree-level value of 0.125. Improved chiral properties 
are seen for the FLIC and MFIC actions, with FLIC4 performing better than 
FLIC12. 

The behavior of the p, nucleon and A masses as a function of squared pion 
mass are shown in Fig. 3 for the various actions. The first feature to note is 
the excellent agreement between the FLIC4 and FLIC12 actions. On the other 
hand, the Wilson action appears to lie somewhat low in comparison. It is also 
reassuring that all actions give the correct mass ordering in the spectrum. The 
value of the squared pion mass at nin/mp = 0.7 is plotted on the abscissa for 
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Table 2. Values of k and the corresponding tt, p, N and A masses for the FLIC action 
with 4 sweeps of smearing at a = 0.7 on a 16^ x 32 lattice at /3 = 4.60. The value for 
KcT is provided in Table 3. A string tension analysis incorporating the lattice coulomb 
term provides a = 0.122(2) fm for = 440 MeV. 



K 


rrin a 


rup a 


mN a 


mA a 


0.1260 
0.1266 
0.1273 
0.1279 
0.1286 


0.5797(23) 
0.5331(24) 
0.4744(27) 
0.4185(30) 
0.3429(37) 


0.7278(39) 
0.6951(45) 
0.6565(54) 
0.6229(65) 
0.5843(97) 


1.0995(58) 
1.0419(64) 
0.9709(72) 
0.9055(82) 
0.8220(102) 


1.1869(104) 
1.1387(121) 
1.0816(152) 
1.0310(194) 
0.9703(286) 



Table 3. Values of k and Ka for the four different actions on a 16^ x 32 lattice at 
/3 = 4.60 which corresponds to a lattice spacing of 0.122(2) fm set by the string tension. 





Wilson 


FLIC12 


FLIC4 


MFIC 




0.1346 


0.1286 


0.1260 


0.1196 


K2 


0.1353 


0.1292 


0.1266 


0.1201 


Ki 


0.1360 


0.1299 


0.1273 


0.1206 


K4 


0.1367 


0.1305 


0.1279 


0.1211 


Ko 


0.1374 


0.1312 


0.1286 


0.1216 




0.1390 


0.1328 


0.1300 


0.1226 



the three actions as a reference point. This point is chosen in order to allow 

comparison of different results by interpolating them to a common value of 
TTiTr/wp = 0.7, rather than extrapolating them to smaller quark masses, which 
is subject to larger systematic and statistical uncertainties. 

The scaling behaviour of the different actions is illustrated in Fig. 4. The 
present results for the Wilson action agree with those of Ref. [35]. The first 
feature to observe in Fig. 4 is that actions with fat-link irrelevant operators 
perform extremely well. For both the vector meson and the nuclcon, the FLIC 
actions perform significantly better than the mean-field improved clover action. 
It is also clear that the FLIC4 action performs systematically better than the 
FLIC12. This suggests that 12 smearing sweeps removes too much short-distance 
information from the gauge-field configurations. On the other hand, 4 sweeps 
of smearing combined with our ©(a^) improved F^;^ provides excellent results, 
without the fine tuning of Csw in the NP improvement program. 

Notice that for the p meson, a linear extrapolation of previous mean-field 
improved clover results in Fig. 4 passes through our mean-field improved clover 
result at a^a ^ 0.075 which lies systematically low relative to the FLIC actions. 
However, a linear extrapolation does not pass through the continuum limit result, 
thus confirming the presence of significant 0{a) errors in the mean-field improved 
clover fermion action. While there are no NP-improved clover plus improved 
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Fig. 3. Masses of the nucleon, A and p meson versus for the FLIC4, FLIC12 and 
Wilson actions on a 16^ x 32 lattice at ,9 = 4.60 which corresponds to a lattice spacing 
of 0.122(2) fm set by the string tension. 



Table 4. String tensions, (3, volumes and results for the vector meson and nucleon 
masses interpolated to mp/mv = 0.7. The scale for the small (3 = 4.60 lattice estimates 
are taken from the large /3 — 4.60 lattice 



p 


Volume 


configs 










4.38 


16-' X 32 


200 


0.371 


2.378(25) 


3.450(35) 


0.8761 


4.53 


20^ X 40 


94 


0.299 


2.318(18) 


3.408(26) 


0.8859 


4.60 


12^ X 24 


200 


0.274 


2.434(26) 


3.554(33) 


0.8889 


4.60 


16^ X 32 


200 


0.274 


2.336(22) 


3.400(26) 


0.8889 


4.80 


16^ X 32 


200 


0.210 


2.427(23) 


3.538(61) 


0.8966 



glue simulation results at a^cr ^ 0.075, the simulation results that are available 
indicate that the fat-link results also compete well with those obtained with a 
NP-improved clover fermion action. 

Having determined FLIC4 is the preferred action, we have increased the 
number of configurations to 200 for this action. As expected, the error bars are 
halved and the central values for the FLIC4 points move to the upper end of 
the error bars on the 50 configuration result, further supporting the promise of 
excellent scaling. 
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Fig. 4. Nucleoli and vector meson masses for the Wilson, NP-improvcd, mean-field 
clover and FLIC actions. Results from the present simulations, based on 50 config- 
urations and indicated by the solid points, are obtained by interpolating the results 
of Fig. 3 to m-nimp — 0.7. The fat links are constructed with n = 4 (solid squares) 
and n — 12 (stars) smearing sweeps at a = 0.7. The solid triangles are results for 
the FLIC4 action when 200 configurations are used in the analysis. The FLIC results 
are offset from the central value for clarity. Our MF clover result at a^a ~ 0.075 lies 
systematically low relative to the FLIC actions. 



In order to further test the scahng of the FLIC action at different lattice 
spacings, we consider four different lattice spacings and three different volumes. 
String tensions, volumes and hadron masses are given in Table 4 and the results 
are displayed in Fig. 5. The two different volumes used at a^a ~ 0.075 indicate 
a small finite volume effect, which increases the mass for the smaller volume at 
a^a ~ 0.075 and ~ 0.045. Examination of points from the small and large vol- 
umes separately indicates continued scaling toward the continuum limit. While 
the finite volume effect will produce a different continuum limit value, the slope 
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Fig. 5. Nucleoli and vector meson masses for the Wilson, 
clover, NP-improved clover and FLIC actions obtained 



Mean-Field (MF) improved 
by interpolating simulation 



results to rriTt/mp = 0.7. For the FLIC action ("FLIC4"), fat links are constructed 
with n = 4 APE-smearing sweeps with smearing fraction a = 0.7, except for the point 
at a^a ~ 0.09 which has n — 6. Results from the current simulations are indicated 
by the solid symbols; those from earlier simulations by open or hatched symbols. The 
solid- lines illustrate fits, constrained to have a common continuum limit, to FLIC, NP- 
improved clover and Wilson fermion action results obtained on physically large lattice 
volumes. 



of the points from the smaller and larger volumes agree, consistent with errors 
oiO{a?). 

Focusing on simulation results from physical volumes with extents ^ 2 fm 
and larger, we perform a simultaneous fit of the FLIC, NP-improved clover 
and Wilson fermion action results. The fits are constrained to have a common 
continuum limit and assume errors arc 0{a?) for FLIC and NP-improved clover 
actions and 0{a) for the Wilson action. An acceptable per degree of freedom 
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Fig. 6. Average number of stabilised biconjugate gradient iterations for the Wilson, 
FLIC and mean-field improved clover (MFIC) actions plotted against m-K/mp. The fat 
links are constructed with n = 4 (solid squares) and n = 12 (stars) smearing sweeps 
at a = 0.7 on a 16^ x 32 lattice at /3 = 4.60 which corresponds to a lattice spacing of 
0.122(2) fm set by the string tension. 



is obtained for both the nucleon and p- meson fits. These results indicate that 
FLIC fermions provide a new form of nonperturbative 0{a) improvement. The 
FLIC fermion results display nearly perfect scaling indicating 0{a^) errors are 
small for this action. 



6 Search For Exceptional Configurations 

Chiral symmetry breaking in the Wilson action allows continuum zero modes of 
the Dirac operator to be shifted into the negative mass region. This problem is 
accentuated as the gauge fields become rough (a large) . Local lattice artifacts 
at the scale of the cutoff (often referred to as dislocations) give rise to spurious 
near zero modes. The quark propagator can then encounter singular behaviour 
as the quark mass becomes light. 

Exceptional configurations are a severe problem in quenched QCD (QQCD) 
because instantons are low action field configurations which appear readily in 
QQCD. These instanton configurations give rise to approximate zero modes 
which should be suppressed at light quark masses by detM which is present 
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Fig. 7. The standard deviation in the error of the n mass for eight quark masses (with 
the star symbols being the lightest quark mass) calculated on 30 configurations plotted 
against the starting configuration number for the FLIC-fermion action on a 20^ x 40 
lattice with a = 0.134(2) fm. 



in the link updates in full QCD. This determinant is not present in QQCD and 
as a result, near-zero modes arc overestimated in the ensemble. 

The addition of the clover term to the fcrmion action broadens the distribu- 
tion of near-zero modes. As a result, the clover action is notorious for revealing 
the exceptional configuration problem in QQCD. The FLIC action is expected to 
reduce the number of exceptional configurations by smoothing the gauge fields of 
the irrelevant operator via APE smearing [5,6]. The smoothing procedure has the 
effect of suppressing the local lattice artifacts and narrowing the distribution of 
near-zero modes, enabling simulations to be performed at light quark masses not 
currently accessible with the standard mean-field or non-perturbative improved 
clover fermion actions. 

In order to access the light quark regime, we would like our preferred action to 
be efficient when inverting the fermion matrix. Fig. 6 compares the convergence 
rates of the different actions on a 16^ x 32 lattice at /3 = 4.60 by plotting the 
number of stabilised biconjugate gradient [34] iterations required to invert the 
fermion matrix as a function of rriTr/mp. For any particular value of TOTr/wp, 
the FLIC actions converge faster than both the Wilson and mean- field improved 
clover fermion actions. Also, the slopes of the FLIC lines are smaller in magnitude 
than those for Wilson and mean-field improved clover actions, which provides 
great promise for performing cost effective simulations at quark masses closer 
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Fig. 8. The n mass calculated for eight quark masses (with the star symbols being the 
lightest quark mass) on 30 configurations plotted against the starting configuration 
number for the FLIC-fermion action on a 20^ x 40 lattice with a = 0.134 fm. 



to the physical values. Problems with exceptional configurations have prevented 
such simulations in the past. 

The ease with which one can invert the fermion matrix using FLIC fermions 
(also see Ref [36]) leads us to attempt simulations down to light quark masses 
corresponding to rriTr/mp = 0.35. Previous attempts with Wilson-style fermion 
actions on configurations with lattice spacing > 0.1 fm have only succeeded in 
getting down to mT^jmp = 0.47 [37]. In order to search for exceptional configura- 
tions, we follow the technique used by Delia Morte et al. [37] and note that in the 
absence of exceptional configurations, the standard deviation of an observable 
will be independent of the number of configurations considered in the average. 
Exceptional configurations reveal themselves by introducing a significant jump 
in the standard deviation as the configuration is introduced into the average. 
In severe cases, exceptional configurations can lead to divergences in correlation 
functions or prevent the matrix inversion process from converging. 

The simulations are on a 20^ x 40 lattice with a lattice spacing of 0.134(2) fm 
set by a string tension analysis incorporating the lattice coulomb term. The 
physical length of the lattice is 2.7 fm. We have used an initial set of 100 
configurations, using n = 6 sweeps of APE-smearing and a five-loop improved 
lattice field-strength tensor. Fig. 7 shows the standard deviation of the pion 
mass for eight quark masses on subsets of 30 (consecutive) configurations with 
a cyclic property enforced from configuration 100 to configuration 1. At first 
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Fig. 9. The standard deviation in the error of the n mass for eight quark masses (with 
the star symbols being the lightest quark mass) calculated on 30 configurations plotted 
against the starting configuration number for the FLIC-fermion action on a 20^ X 40 
lattice with a = 0.134 fm. Configuration numbers 2, 13, 30, 34, 41 and 53 have been 
omitted. 



glance, it is obvious that the error blows up for several quark masses at A'' = 
12 and drops again at A'' = 42. As configurations 12 through 41 are included 
in the average at A'^ = 12, this indicates that configuration number 41 is a 
candidate for an exceptional configuration. An inspection of the pion mass in 
Fig. 8 shows that the pion mass for the third lightest quark mass decreases 
significantly more than the second or fourth lightest quark masses. This indicates 
that Kcr for this configuration lies somewhere between kq and K7. A solution to 
this problem would be to use the the modified quenched approximation (MQM) 
from Ref. [22] and move on this configuration back to the ensemble average for 
KcT- However, since the movement of Kct is largely a quenched artifact and would 
be suppressed in a full QCD simulation we prefer to simply identify and remove 
such configurations from the ensemble. Obviously, if we find that a significant 
percentage of our configurations are having trouble at a particular quark mass, 
then it would make no sense to proceed with the simulation. We would then 
have to conclude that we have reached the light quark mass limit of our action 
and simply step back to the next lightest mass. 

Now let us return to Fig. 7. In addition to the highly exceptional configu- 
ration number 41, we also notice a large increase in error in the lightest quark 
mass for configuration numbers 2, 13, 30, 34 and 53. Upon removal of these con- 
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Table 5. Values of k and the corresponding tt, p, N and A masses on a 20^ x 40 lattice 
for the FLIC action with 6 sweeps of smearing at ot = 0.7. A string tension analysis 
incorporating the lattice coulomb term provides a = 0.134(2) fm for y/a = 440 MeV. 





a 


rUp a 


m-K a 


rriA a 


0.1278 


0.5400(30) 


0.7304(55) 


1.0971(80) 


1.2238(98) 


0.1283 


0.4998(31) 


0.7053(58) 


1.0522(84) 


1.1899(102) 


0.12885 


0.4521(34) 


0.6774(63) 


1,0006(91) 


1.1528(108) 


0.1294 


0.3990(38) 


0.6491(72) 


0.9465(101) 


1.1162(115) 


0.1299 


0.3434(43) 


0.6228(87) 


0.8944(116) 


1.0841(125) 


0.13025 


0.2978(47) 


0.6040(107) 


0.8562(134) 


1.0630(135) 


0.1306 


0.2419(54) 


0.5845(143) 


0.8172(171) 


1.0443(154) 


0.1308 


0.1972(69) 


0.5812(213) 


0.7950(215) 


1.0380(189) 



figurations, we see in Fig. 9 a near-constant behaviour of the standard deviation 
for the remaining configurations. This means that our ehmination rate for our 
FLIC6 action on a lattice with a spacing 0.134 fm is about 6%. So for the 100 
configurations used in this analysis, we are able to use 94 of them to extract 
hadron masses. 

A similar analysis on a 16^ x 32 lattice at /? = 4.60 providing a finer lattice 
spacing of 0.122(2) fm reveals a much smaller exceptional configuration rate. In 
a sample of 200 configurations, 4 were identified as exceptional. The increase 
from 2% to 6% in going from a ~ 0.125 to 0.135 fm suggests that the coarser 
lattice spacing is near the limit of applicability for FLIC fermions in the light 
quark mass regime. 

7 Octet-Decuplet Mass Splittings 

The results presented in this section are based on an initial sample of 94 gauge- 
field configurations of an anticipated 400 configurations. Figure 10 shows the 
N and A masses as a function of for the FLIC-fermion action on 20^ x 40 
lattices with a = 0.132 fni (which corresponds to a string tension scale with 
^ = 450 MeV) such that the nucleon extrapolation passes through the physical 
value for clarity. An upward curvature in the A mass for decreasing quark mass is 
observed in the FLIC fermion results. This behaviour, increasing the quenched 
N — A mass spitting, was predicted by Young et al. [38,39] using quenched 
chiral perturbation theory (QxPT) formulated with a finite-range regulator. A 
fit to the FLIC-fermion results is illustrated by the solid curves. The dashed 
curves estimate the correction that will arise in unquenching the lattice QCD 
simulations [38,39]. We note that after we have corrected for the absence of sea 
quark loops, our results agree simultaneously with the physical values for both 
the nucleon and A. 

We also calculate the light quark mass behaviour of the octet and decuplet 
hyperons. The strange quark mass is chosen in order to reproduce the physical 




Fig. 10. Nucleon and A masses for the FLIC-fermion action on a 20^ x 40 lattice. Here 
we select a = 0.132 fm (which corresponds to the string tension with = 450 MeV) 
such that the nucleon extrapolation passes through the physical value for clarity. The 
solid curves illustrate fits of finite-range regularised quenched chiral perturbation theory 
[38,39] to the lattice QCD results. The dashed curves estimate the correction that will 
arise in unquenching the lattice QCD simulations [38,39]. Stars at the physical pion 
mass denote experimentally measured values. 



strange quark mass according to the phenomenological value of an ss pseu- 
doscalar meson, 

mis = 2ml -ml. (41) 

Upon substitution of the physical masses for the it and K mesons, this corre- 
sponds to an ss pseudoscalar meson mass of ~ 0.470 GeV^ which occurs at our 
third heaviest quark mass in the 20^ x 40 lattice analysis. The results from this 
calculation are given in Table 6 and are illustrated in Fig. 11. The results show 
the correct ordering and in particular, we notice a mass splitting between the 
strangeness = —1 (I = I) E and {I = 0) A baryons becoming evident in the 
light quark mass regime. 

Just as we saw the; non-analytic behaviour of quenched chiral perturba- 
tion theory in the Zi-baryon mass in Fig. 10 leading to an enhancement of 
the quenched N — A mass spitting. Fig. 12 shows a similar enhancement for 
the decuplet-octet mass splittings in S and S baryons respectively. The quark 
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Fig. 11. Octet (top) and decuplet (bottom) baryon masses for the FLIC-fermion action 
on a 20^ x 40 lattice witli a = 0.134 fm. 
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Fig. 12. Octet and decuplet baryon masses for S (top) and ^ (bottom) for the FLIC- 
fermion action on a 20^ x 40 lattice with a = 0.134 fm. 
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Table 6. Values of k, the octet A, E, S and decuplet E* , S* masses on a 20^ X 40 
lattice for the FLIC action with 6 sweeps of smearing at a = 0.7. A string tension 
analysis provides a = 0.134(2) fm for ^/a = 440 MeV. 



K 


nijx a 


nis a, 


rn= a 


rris , a 


rriS f a 


0.1278 


1.0696(84) 


1.0616(83) 


1.0381(87) 


1.2002(101) 


1.1765(104) 


0.1283 


1.0376(86) 


1.0328(86) 


1.0206(88) 


1.1776(104) 


1.1652(106) 


0.12885 


1.0006(91) 


1.0006(91) 


1,0006(91) 


1.1528(108) 


1.1528(108) 


0.1294 


0.9615(97) 


0.9680(97) 


0.9799(94) 


1.1284(113) 


1.1406(110) 


0.1299 


0.9235(106) 


0.9383(106) 


0.9603(98) 


1.1070(118) 


1.1299(113) 


0.13025 


0.8955(117) 


0.9178(116) 


0.9462(102) 


1.0930(124) 


1.1229(115) 


0.1306 


0.8667(137) 


0.8980(132) 


0.9323(109) 


1.0806(132) 


1.1166(118) 


0.1308 


0.8544(154) 


0.8919(152) 


0.9254(114) 


1.0772(142) 


1.1145(120) 



model predicts that the hypcrfine splittings should approximately satisfy [40] 

S:-E,=n,tXg = S:-E,, (42) 

where the baryon label denotes the hypcron mass and fig (fiq) denotes the mag- 
netic moment of the strange (light) constituent quark. Fig. 13 shows that even 
though the quenched approximation enhances the splitting between octet and 
decuplet baryons, the splittings for the S and S baryons still satisfy Eq. (42). 
Agreement of the quenched QCD results with the quark model prediction is 
not surprising since both have a suppressed meson cloud. Similarly, one expects 
further suppression of the meson cloud when two (heavy) strange quarks are 
present in a baryon. 

8 Summary 

We have calculated hadron masses to test the scaling of the Fat-Link Irrelevant 
Clover (FLIC) fermion action, in which only the irrelevant, higher-dimension 
operators involve smeared links. One of the main conclusions of this work is 
that the use of fat links in the irrelevant operators provides a new form of 
nonperturbative 0{a) improvement. This technique competes well with 0{a) 
nonpcrturbativc improvement on mean field-improved gluon configurations, with 
the advantage of a reduced exceptional configuration problem. 

Quenched simulations at quark masses down to m^/mp = 0.35 have been 
successfully performed on a 20"^ x 40 lattice with a lattice spacing of 0.134(2) fm 
on 94 out of 100 configurations. Simulations at such light quark masses reveal 
the non-analytic behaviour of quenched chiral perturbation theory and provide 
for an interesting analysis of the hyperfine splittings between octet and decuplet 
baryons. 
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Fig. 13. Docuplct [Md) - octet (Mo) baryon mass splittings for the FLIC-fermion 
action on a 20^ x 40 lattice with a = 0.132 fm. 
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